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ABSTRACT 

We present high-speed ULTRACAM and SOFI photometry and X-shooter spec- 
troscopy of the recently discovered 94 minute orbital period eclipsing white dwarf / 
main-sequence binary SDSS J085746. 18+034255. 3 (CSS 03170) and use these observa- 
tions to measure the system parameters. We detect a shallow secondary eclipse and 
hence are able to determine an orbital inclination of i = 85.5° ± 0.2°. The white dwarf 
has a mass of 0.51 ± O.O5M and a radius of 0.0247 ± 0.0008Rq. With a temperature 
of 35, 300 ± 400K the white dwarf is highly over-inflated if it has a carbon-oxygen core, 
however if it has a helium core then its mass and radius are consistent with evolution- 
ary models. Therefore, the white dwarf in SDSS J085746. 18+034255.3 is most likely a 
helium core white dwarf with a mass close to the upper limit expected from evolution. 
The main-sequence star is an M8 dwarf with a mass of 0.09 ± O.OIMq and a radius of 
0.110 ± 0.004Rq placing it close to the hydrogen burning limit. The system emerged 
from a common envelope ~ 20 million years ago and will reach a semi-detached con- 
figuration in ~ 400 million years, becoming a cataclysmic variable with a period of 66 
minutes, below the period minimum. 

Key words: binaries: eclipsing - stars: fundamental parameters - stars: low mass - 
white dwarfs - stars: individual: SDSS J085746. 18+034255.3 



INTRODUCTION and radii are not possible. Despite their importance, the 

mass-radius relations for both white dwarfs and low-mass 



stars remain largely untested. 



Mass-radius relations are of fundamental importance in a 
wide range of astrophysical circumstances. They are rou- 
tinely used to infer accurate masses and radii of transit- IProvencal et~aH (|l998l ) tested the white dwarf mass- 
ing exoplanets, calibrating stellar evolutionary models and radius relation using Hipparcos parallaxes to determine the 
understanding the late evolution of mass transferring bi- radii for white dwarfs in visual binaries, common proper- 
naries such as cata clysmic variables (|Littlefair et all |200S|; motion (CPM) systems and field white dwarfs. However, 
ISavoury et aD|201l|). Additionally, the mass-radius relation the radius measurements for all of these systems still rely 
for white dwarfs has played an im portant role in estim ating to some extent on model atmosphere calculations. For field 
the distance to globular clusters dRenzini et al]|l99 6|) and white dwarfs t he m ass determinations are also indirect. 
the d etermination of the age of the Galactic disk (|Wood| iBarstow et ail (|2005l ) used Hubble Space Telescope/STIS 
1 1992!) . Mass-radius relations play a crucial role in determin- spectra to measure the mass of Sirius B to high precision, 
ing stellar parameters of single, isolated stars as well as in however, their radius constraint still relied on model atmo- 
non-eclipsing binaries, where direct measurements of masses sphere calculations and is therefore less direct when it comes 

to testing white dwarf mass-radius relations. To date, only 
two white dwarfs have had their mas ses and radii model - 

* steven.parsons@warwick.ac.uk independently measured, V471 Tau (|0'Brien et al.l 1200 if ) 
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Table 1. Journal of observations. Exposure times for X-shooter observations are for UVB arm, VIS arm and NIR arm respectively. The 
primary eclipse occurs at phase 1, 2 etc. 



Date at 


Instrument 


Filter (s) 


Start 


Orbital 


Exposure 


Conditions 


start of run 






(UT) 


phase 


time (s) 


(Transparency, seeing) 


2010/12/03 


ULTRACAM 


u'g'r' 


07:14 


0.86-1.08 


4.0 


Excellent, ~ 1 arcsec 


2010/12/07 


ULTRACAM 


u'g'r' 


06:53 


0.95-1.20 


4.0 


Excellent, ~ 1 arcsec 






u'g'r' 


UU.Ol 


n 9i n fi^ 

U.Zl u.oo 


r , n 

■ ) . u 


Excellent, r ^ 1 arcsec 


2010/12/09 


ULTRACAM 


u'g'i' 


06:49 


0.77-1.80 


5.0 


Excellent, ~ 1 arcsec 


2010/12/10 


ULTRACAM 


u'g'r' 


06:31 


0.93-1.29 


4.0 


Good, 1.0 — 1.5 arcsec 


2010/12/16 


ULTRACAM 


u'g'i' 


04:16 


0.68-2.12 


4.0 


Good, 1.0 — 2.0 arcsec 


2011/01/07 


ULTRACAM 


u'g'i' 


06:53 


0.54-1.11 


4.0 


Good, 1.0 — 1.5 arcsec 


2011/01/09 


ULTRACAM 


u'g'r' 


04:46 


0.70-1.69 


4.0 


Good, 1.0 — 1.5 arcsec 


2011/01/10 


ULTRACAM 


u'g'r' 


04:42 


0.03-0.85 


4.0 


Good, 1.0 — 1.5 arcsec 


2011/01/11 


ULTRACAM 


u'g'r' 


03:22 


0.53-1.95 


4.0 


Good, 1.0 — 1.5 arcsec 


2011/02/27 


X-shootcr 




03:59 


0.94-2.14 


191,213,480 


Good, 1.0 — 1.5 arcsec 


2011/03/03 


X-shooter 




02:15 


0.27-1.35 


191,213,480 


Good, 1.0 — 1.5 arcsec 


2011/04/06 


SOFI 


J 


00:58 


0.72-1.27 


15.0 


Excellent, < 1 arcsec 



and NN Ser (|Parsons et al.ll2~010al ). Both of these systems are 
eclipsing post common envelope binaries (PCEBs), demon- 
strating the importance of eclipsing systems in determining 
precise masses and radii. 

PCEBs are the product of main sequence binary stars 
where the more massive members of the binaries evolve off 
the main-sequence and fill their Roche lobes on either the 
giant or asymptotic giant branches. This initiates a dynami- 
cally unstable mass transfer to the less massive components 
which are unable to accrete the material, hence common en- 
velopes of material are formed containing the core of the 
giant and the companion star. Frictional forces within this 
envelope cause the two stars to spiral inwards. The ensu- 
ing loss of angular momentum expels the common envelope 
revealing the tightly bound core and companion. An addi- 
tional advantage to studying PCEBs is that some contain 
low-mass helium core white dwarfs, enabling us to measure 
masses and radii for these white dwarfs which are not found 
in wide binaries (that did not undergo a common envelope 
phase) . 

PCEBs usually contain a low-mass late- type com panion 
to the white dwarf (|Rebassa-M anscrg as et al.ll201ll ). Hence 
an additional benefit of studying PCEBs is that, under 
favourable circumstances, not only are the white dwarf's 
mass and radius determined independently of any model, 
so too are the mass and radius of its companion. There is 
disagreement between models and observations of low mass 
stars; the models t end to under predict the radii by as much 
as 20-30 per cent (|L6pez- Mor ales 2007). However, it is un- 
clear whether this effect is a genuine problem in the models, 
or a systematic effect caused by testing the mass-radius re- 
lation with close binaries, which are tidally locked into rapid 
rotation and thus can be inflated by stronger ma gnetic fields 
inhibi ting the efficiency of convection. Recently. iKraus et al.l 
(2011) analysed 6 new M dwarf eclipsing binary systems and 
found that the stars in short period binaries (P < 1 day) 
have radii that are inflated by up to 10 per cent compared 
to evolutionary models, whilst those in longer period bina- 
ries have smaller radii. This supports the inflation via rapid 
rotation hypothesis. However, their sample contains stars of 
masses M > OAMq and thus does not test this effect in the 
very low mass range. Indeed the radius of the rapidly rotat- 



ing low- mass (M = O.111M ) M dwarf in the PCEB NN 
Ser showed very little inflation above evolutionary models, 
and what inflation there was could be entirely explained as 
the result of irradiation from its hot w hite dwarf compan - 
ion (| Parsons et a l. 2010a). Furthermore. lBerger et al] (2006) 
found a similar underestimation of the radius of single stars 
compared to evolutionary models, implying that rapid rota- 
tion may not be the only cause, although they note that 
their lower metallicity targets do show better agreement 
with models. 

The subject of this paper, SDSS J085746. 18+034255.3 
(henceforth SDSS J0857+0342), was first discovered as a hot 
(Tea = 36181 ± 390K) white dwarf from the Sloan Digital 
Sky Survey byjEiscnstcin et al] (|2006l ). Regular eclipses were 
observed bv lDrake et al.l (|2010l ) as part of the Catalina Real- 
time Transient Survey. They determined an orbital period 
of 94 minutes making SDSS J0857+0342 the shortest known 
period, detached, white dwarf / main-sequence binary. 

In this paper we present a combination of ULTRA- 
CAM and SOFI high-speed photometry and X-shooter spec- 
troscopy of SDSS J0857+0342 and use these to constrain the 
system parameters. We also discuss the past and future evo- 
lution of the system. 

2 OBSERVATIONS AND THEIR REDUCTION 

2.1 NTT/ULTRACAM photometry 

SDSS J0857+0342 was observed with ULTRACAM 
mounted as a visitor instrument on the New Technology 
Telescope (NTT) at La Silla between December 2010 and 
January 2011. ULTRACAM i s a high-speed, triple-beam 
CCD camera l|Dhillon et al] 120071 1 which can acquire 
simultaneous images in three different bands; for our obser- 
vations we used the SDSS u , g' and either r' or il filters. A 
complete log of these observations is given in Table \T\ We 
windowed the CCD in order to achieve exposure times of 
~4 seconds which we varied to account for the conditions; 
the dead time between exposures was ~ 25 ms. 

All of these data were reduced using the ULTRACAM 
pipeline software. Debiassing, flatfielding and sky back- 
ground subtraction were performed in the standard way. 
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Figure 1. Spectral model fit to the white dwarf component of SDSS J0857+0342 from the X-shooter spectroscopy. Left: best fit (black 
lines) to the normalised H/3 to H9 line profiles (gray lines, top to bottom). Interstellar calcium absorption is also visible. Right: 1, 2 and 
3cr x 2 contour plots in the T e g — logp plane to the line profile fits. The intrinsic contribution from the secondary star in this wavelength 
range is negligible. We remove the reflection component by only selecting those spectra taken around the primary eclipse, where the 
reflection is minimal. 



The source flux was determined with aperture photome- 
try using a variable aperture, whereby the radius of the 
aperture is scaled according to the full width at half maxi- 
mum (FWHM). Variations in observing conditions were ac- 
counted for by determining the flux relative to a comparison 
star in the field of view. The data were flux calibrated by 
determining atmospheric extinction coefficients in each of 
the bands in which we observed and we calculated the abso- 
lute fl ux of our target u sing observations of standard stars 
(from lSmith et aTll2002l ) taken in twilight. The comparison 
star used (at 08:57:59.04 +03:41:54.5) has magnitudes of 
u = 16.209, g' = 14.997, r' = 14.602 and i' = 14.455. Using 
our extinction coefficients we extrapolated all fluxes to an 
airmass of 0. The systematic error introduced by our flux 
calibration is < 0.1 mag in all bands. 

2.2 VLT/X-shooter spectroscopy 

We obtained servi ce mode observat i ons o f SDSS J0857+0342 
with X-shooter jD'Odorico et al l 120061 ) mounted at the 
VLT-UT2 telescope. The observations were designed to 
cover two entire orbits of the system. Details of these ob- 
servations are listed in Table [1] X-shooter is a medium res- 
olution spectrograph consisting of 3 independent arms that 
give simultaneous spectra longward of the atmospheric cut- 
off (0.3 microns) in the UV (the "UVB" arm), optical (the 
"VIS" arm) and up to 2.5 microns in the near-infrared (the 
"NIR"arm). We used slit widths of 0.8", 0.9" and 0.9" in 
X-shooter's three arms and 2x2 binning resulting in a res- 
olution of R ~ 6, 000. After each exposure we nodded along 
the slit to help the sky subtraction in the NIR arm. 

The reduction of the raw frames was conducted us- 
ing the standard pipeline release of the X-shooter Common 
Pipeline Library (CPL) recipes (version 1.3.0) within ES- 
ORex, the ESO Recipe Execution Tool, version 3.9.0. The 
standard recipes were used to optimally extract and wave- 
length calibrate each spectrum. The instrumental response 
was removed by observing the spectrophotometric standard 
star EG 274 and dividing it by a flux table of the same star 



to produce the response function. We then heliocentrically 
corrected the wavelength scales of each of the spectra. For 
the NIR arm we combine frames taken at different nod po- 
sitions to improve the sky subtraction. However, given the 
long exposure times in the NIR arm, this leads to a dramatic 
loss of orbital phase resolution. We achieve a signal-to-noise 
(S/N) of ~ 20 in the UVB arm per exposure, ~ 10 in the 
VIS arm per exposure and ~ 10 in the NIR in two hours. 
The K band data are completely dominated by sky noise. 
The S/N in the NIR arm is too low to detect any features 
from the secondary star. 

The ULTRACAM and SOFI light curves were used to 
flux calibrate the X-shooter spectra. We fitted a model to 
each of the ULTRACAM light curves (see Section 13. 2|) and 
the SOFI J-band light curve in order to reproduce the light 
curve as closely as possible. The model was then used to 
predict the flux at the times of each of the X-shooter obser- 
vations. We then derived synthetic fluxes from the spectra 
for the ULTRACAM u',g',r', i! and z filters as well as the 
SOFI J and H filters. We extrapolated the light curve mod- 
els to those bands not covered by our photometry. We then 
calculated the difference between the model and synthetic 
fluxes and fitted a second-order polynomial to them. This 
correction was then applied to each spectrum. This corrects 
for variable extinction across the wavelength range, as well 
as variations in seeing. 

2.3 NTT/SOFI J band photometry 

SDSS J0857+0342 was observed with SOFI mounted at the 
NTT on the 6th April 2011. The observations were made 
in fast photometry mode equipped with a J band filter. We 
windowed the detector to achieve a cycle time of 15 seconds 
and offset the telescope every 10 minutes in order to improve 
sky subtraction. 

Debiassing (which also removes the dark current) and 
flatfielding were performed in the standard way. Sky sub- 
traction was achieved by using observations of the sky when 
the target had been offset. The average sky level was then 
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Figure 2. Top to bottom: phase-folded ULTRACAM u'g'r'i' and SOFI J band light curves of SDSS J0857+0342. The ULTRACAM 
light curves have been binned by a factor of three. Each light curve is offset by 0.2 mjy from the one below it. The J band light curve 
has been scaled up by a factor of three for clarity. 



added back so that we could determine the source flux and 
its uncertainty with standard aperture photometry, using a 
variable aperture, within the ULTRACAM pipeline. A com- 
parison star (2MASS 08575095+0340301, J = 13.289) was 
used to account for variations in observing conditions. Flux 
calibration was done using the comparison s tar J band mag- 
nitud e retrieved from the 2MASS catalogue (Skrutskie et al.l 

l2ooeh . 



3 RESULTS 

3.1 White dwarf temperature 

lEisenstein et al.1 \200dh and brake et all |2010f i both 
determined the temperature of the white dwarf in 
SDSS J0857+0342 using its SDSS spectra. However, the re- 
processed light from the heated face of the secondary star 
can lead to incorrect results from spectral fitting. There- 
fore, we combined the X-shooter spectra taken around the 
eclipse, where the reflection effect is minimal, and used 
this spectrum to determine the t emperature and log. q val- 
ues us ing the method outlined in iRebassa-Mansergas et al.l 
(2007). Figure [T] shows the results of the spectral fit to 
the white dwarf features, we determined a temperature of 
35,300 ± 400K and a logg of 7.3 ± 0.1, which is consis- 
tent with the ligh t curve model f i ts an d th e fits to the 
SPSS spectra from lEisenstein et al.l (|2006t ) and lDrake et al.l 
|2010f ). The fit also gave a distance to SDSS J0857+0342 of 
968 ± 55pc, consistent with previous measurements. 



Table 2. NTT/ULTRACAM mid-eclipse times for 
SDSS J0857+0342. 



Cycle 


Eclipse time 


Uncert 


Number 


MJD(BTDB) 


MJD(BTDB) 


-298 


55533.3139971 


0.0000024 


-237 


55537.2848843 


0.0000019 


-206 


55539.3028807 


0.0000024 


-191 


55540.2793255 


0.0000022 


-100 


55546.2031193 


0.0000064 


-99 


55546.2682076 


0.0000020 


240 


55568.3359334 


0.0000020 


269 


55570.2237344 


0.0000024 


299 


55572.1766308 


0.0000025 



3.2 Light curve model fitting 



for 



Figure [2] shows phase folded light 
SDSS J0857+0342 in the ULTRACAM ugri 
and the SOFI J band. The deep primary eclipse is clearly 
visible in all bands; flux is detected in eclipse during the J 
band eclipse only. A strong reflection effect is also evident 
in each of the bands caused by reprocessed light from the 
irradiated face of the secondary star. Also visible is the 
shallower secondary eclipse which becomes deeper in the 
redder bands (this phase was not covered by the SOFI 
observations). The detection of the secondary eclipse allows 
us to determine the inclination of the system and establish 
precise scaled radii, R sec /a and Rwo/a, where a is the 
orbital separation. 

To measure the system parameters we fitted the UL- 
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Orbital Phase Orbital Phase 

Figure 3. Model fits to the ULTRACAM light curves with residuals shown below. Top: Full orbit in the i' band. Bottom left: The g' 
band primary eclipse. Bottom right: The r' band secondary eclipse, also shown is the same model but without the secondary eclipse. 



TRACAM light curves using a code written to produce mod- 
els f or the general case of bin aries containing white dwarfs 
fsee ICopperwheat et al.l 120101 for details). It has been used 
in the study of other white dwarf-main sequence binaries 
l|Pvrzas et al.ll2009l ; [Parsons et al.l2010bl l. The program sub- 
divides each star into small elements with a geometry fixed 
by its radius as measured along the direction of centres to- 
wards the other star. Roche geometry distortion and irra- 
diation of the secondary star are included, the irradiation 
is approximated by o-T'* cc — tjT s * c + AFi TT where T' sec IS 
the modified temperature and T scc is the temperature of the 
unirradiated main-sequence star, a is the Stefan-Boltzmann 
constant, A is fraction of the irradiating flux from the white 
dwarf absorbed by the secondary star and Fi lr is the irradi- 
ating flux, accounting for the angle of incidence and distance 
from the white dwarf. 

We used the Markov Chain Monte Carlo (MCMC) 
method t o determine the distributions of our model pa- 
rameters (|Press et all 120071 ). The MCMC method involves 
making random jumps in the model parameters, with new 
models being accepted or rejected according to their proba- 
bility computed as a Bayesian posterior probability. In this 
instance this probability is driven by x 2 with all prior prob- 
abilities assumed to be uniform. A crucial practical consid- 
eration of MCMC is the number of steps required to fairly 
sample the parameter space, which is largely determined by 



how closely the distribution of parameter jumps matches the 
true distribution. We therefore built up an estimate of the 
correct distribution starting from uncorrelated jumps in the 
parameters, after which we computed the covariance matrix 
from the resultant chain of parameter values. The covari- 
ance matrix was then used to define a multivariate normal 
distribution that was used to make the jumps for the next 
chain. At each stage the actual size of the jumps was scaled 
by a single factor set to deliver a model acceptance rate 
of ~ 25 per cent. After 2 such cycles, the covariance ma- 
trix showed only small changes, and at this point we carried 
out the long "production runs" during which the covariance 
and scale factor which define the parameter jumps were held 
fixed. Note that the distribution used for jumping the model 
parameters does not affect the final parameter distributions, 
only the time taken to converge towards them. 

The parameters needed to define the model were: the 
mass ratio, q = M sbc /Mwd, the inclination, i, the scaled 
radii R SB c/a and Rwo/a, the unirradiated temperatures, 
Teff,wD and T e s jSe c! limb darkening coefficients for the both 
stars, the time of mid eclipse, Tb and the period, P. 

Each light curve was initially fitted in order to deter- 
mine the mid-eclipse times and thus an accurate ephemeris. 
The mid-eclipse times are listed in Table [2] We update the 
ephemeris to 
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Table 3. Parameters from Markov chain Monte Carlo minimisation, some fitted, some fixed a priori. a^D an d £>wd are the quadratic 
limb darkening coefficients of the white dwarf, a Bec is the linear limb darkening of the secondary star. A is the fraction of the irradiating 
flux from the white dwarf absorbed by the secondary star. 



Parameter 


v! 


g' 


r' 


i' 


i 


86.6 ±0.5 


85.4 ±0.2 


85.5 ±0.2 


84.8 ±0.5 




0.0461 ± 0.0012 


0.0438 ± 0.0006 


0.0431 ± 0.0005 


0.0423 ± 0.0013 


r scc /a 


0.180 ±0.004 


0.192 ±0.003 


0.191 ±0.003 


0.203 ± 0.008 


-^eff , sec 


4922 ± 97 


3575 ± 67 


3298 ± 64 


3016 ± 104 


"WD 


0.089 


0.064 


0.057 


0.054 


bwu 


0.199 


0.178 


0.143 


0.119 


&SCC 


-0.90 ±0.10 


-0.30 ±0.05 


-0.16 ± 0.05 


-0.09 ± 0.11 


A 


0.83 ±0.02 


0.64 ±0.02 


0.78 ±0.02 


0.69 ± 0.06 



MJD(BTDB) = 55552.712 76517(78) + 0.065 096 538(3)£, 



which is consistent with the ephemeris of IDrake et alj 
l|2010l ). 

We then phase-folded the data and kept the period fixed 
but allowed the time of mid eclipse to vary. We kept the mass 
ratio fixed at 0.13 (the light curves are only weakly depen- 
dent on this parameter) and the temperature of the white 
dwarf fixed at 35, 300K (see section f3 . 1 [) . We determined 
the quadratic limb darkening coefficients of the white dwarf 
based on a white dwarf atmosphere model with a tempera- 
ture of T e B — 35, 300K and logg = 7.3 folded through th e 
ULTRACAM u'g'r'i' filter profiles (|Gansicke et al.lll995h . 
The coefficients (owd and 6wd) for each filter are listed in 
Table[3]for J(/i)//(l) = 1 - a(l - y) - 6(1 - fif, where [i is 
the cosine of the angle between the line of sight and the sur- 
face normal; we kept these values fixed. All other parameters 
were allowed to vary and several chains were started with 
different initial parameters in order to check if the results 
were consistent and that we had found the global minimum. 

The best fit parameters and their statistical errors are 
displayed in Table Figure [3] shows the fits to the light 
curves over the entire orbit as well as zoomed in around the 
two eclipses. The MCMC chains showed no variation beyond 
that expected from statistical variance. Since we didn't cover 
the secondary eclipse in the J band and the S/N of the SOFI 
data is quite low, we don't fit this data. 

The four different bands give consistent results, al- 
though the u' band favours a slightly higher inclination 
but is less constrained due to the shallowness of the sec- 
ondary eclipse in this band. The faintness of the system in 
the i! band limits the constraints in this band. We found 
an inclination of 85.5° ± 0.2° and scaled radii of Rwo/a = 
0.0436±0.0004 and R scc /a = 0.190±0.002. Given our black 
body assumption, T sec only approximately represents the 
true temperature of the secondary star, its more important 
role is as a flux scaling factor. An interesting trend is seen 
in the limb darkening coefficients for the secondary star, 
which are all negative (limb brightening) and the amount of 
limb brightening decreases with increasing wavelength. This 
is presumably the result of seeing to different depths at dif- 
ferent wavelengths; a similar effect was seen with NN Ser 
jParsons et al.ll2010al ). 



Table 4. Hydrogen Balmer line offsets and velocities. 



Line 


7WD 


%D 


7cmis 


^cmis 




(kms -1 ) 


(kms -1 ) 


(kms -1 ) 


(kms- 1 ) 


H 6 






28.9 ±4.9 


319.4 ±8.7 


H 7 


48.4 ±8.6 


65.9 ± 11.8 


25.5 ±5.3 


324.7 ±8.9 


H/3 


50.6 ±5.8 


61.6 ±8.8 


28.5 ±5.3 


343.6 ±9.8 


H a 


42.1 ±6.3 


66.7 ± 10.0 


31.8 ±5.0 


336.2 ±9.9 



3.3 Radial Velocities 

Figure [4] shows an average spectrum of SDSS J0857+0342. 
The spectral features seen are similar to other PCEBs, 
namely hydrogen Balmer absorption originating from the 
white dwarf and emission from the heated face of the sec- 
ondary star, which disappears around the primary eclipse. 
Calcium emission is also seen as well as deep interstellar 
calcium H and K absorption which show no radial velocity 
variations; no other features from either star are visible. Fig- 
ure [5] shows trailed spectra of the four longest wavelength 
hydrogen Balmer lines. The emission from the secondary 
star can be seen as well as the non-LTE absorption core 
from the white dwarf moving in anti-phase. 

For each of the Balmer lines we fit both the absorption 
and emission components together by simultaneously fitting 
all of the spectra. We used this approach because the S/N 
on an individual spectrum was quite low and the absorption 
from the white dwarf is filled in by the emission compo- 
nent from the secondary star around phase 0.5. Therefore, 
by fitting all of the spectra simultaneously we obtain a more 
robust estimate of the radial velocities. We used a combi- 
nation of a straight line and Gaussians for each spectrum 
(including a broad Gaussian component to account for the 
wings of the primary star's absorption) and allowed the po- 
sition of the Gaussians to change velocity according to 

V = 7 + ATsin(27T0), 

for each star, where 7 is the velocity offset of the line from its 
rest wavelength and (j> is the orbital phase of the spectrum. 
We allow the height of the emission component to vary with 
orbital phase according to 

H =H - Hi cos(2tt0), 

which allows the height to peak at phase 0.5, where the 
irradiation effect is largest. This approach gave better fits 
than keeping the height at a fixed value. 



© 2011 RAS, MNRAS 000,[THTT] 



The Eclipsing binary SDSS J085746. 18+034255.3 7 




0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 

Wavelength (microns) 



Figure 4. Averaged X-shooter spectrum of SDSS J0857+0342. No features from the secondary star are seen in the NIR, the only features 
visible are telluric or sky lines that have not been completely removed. Those spectra taken during the eclipse were not included in the 
averaging. Inset is a zoom in on the upper Balmer series, interstellar calcium absorption is also visible. 



Table |4] lists the offsets and radial velocities for the four 
longest wavelength hydrogen Balmer lines, the shorter wave- 
length lines were unsuitable for fitting since the emission 
component is much weaker and the absorption components 
lack the non-LTE core, this is also the case for the H S 
absorption. The results for each line are consistent and we 
adopt values of -Kwd = 64 ± 6 kms" 1 and A" cm i s = 330 ± 5 
kms" 1 . We also measure a gravitational redshift for the 
white dwarf of zwd = 7wd — 7cmis = 18 ± 5 kms" 1 . 



3.4 ivscc correction 

The emission lines seen in the X-shooter spectra are the re- 
sult of reprocessed light from the surface of the secondary 
star facing the white dwarf, hence their radial velocity am- 
plitude represents a lower limit to the true centre of mass 
radial velocity amplitude. For accurate mass determinations 
the centre of mass radial velocity amplitude is required, thus 
we need to determine the deviation between the reprocessed 
light centre and the centre of mass for the secondary star. 
The radial velocity of the centre of mass (K scc ) is related to 
that of the emission lines (K eln i a ) by 

K scc = Kcmi \ R (1) 

1-/(1 + 9)^ 

where / is a constant between and 1 which depends upon 
the location of the centre of light. For / = the emission is 
spread uniformly across the entire surface of the secondary 
star and therefore the centre of light is the same as the centre 
of mass. For / = 1 all of the flux is assumed to come from 
the point on the secondary star's surface closest to the white 
dwarf (the substellar point). The maximum possible value 
for K Bec is therefore 423 kms" (using A'cmis = 330 kms" 1 , 
R S cc/a = 0.190 and / = 1). 

The centre of light for an emission line is related to the 
optical depth of the emission (|Parsons et al.ll2010aT ). Opti- 
cally thick emission tends to be preferentially radiated per- 
pendicular to the stellar surface, therefore at the quadra- 
ture phases, we will see the limb of the irradiated region 
more prominently (compared to the region of maximum 
irradiation) than we would otherwise. This will lead to a 



higher observed semi-amplitude and hence a smaller correc- 
tion factor is needed. The reverse is true for optically thin 
lines where the emission is radiated equally in all directions, 
hence emission from the substellar point becomes enhanced 
at quadrature, leading to a low semi-amplitude and a larger 
correction factor. The correction factor for an optically thin 
line (assumi ng emissivity proport ional to the incident flux) 
is / = 0.77 {P arsons et alj|2010ah . Therefore, a more strin- 
gent upper limit on the correction factor can be set by us- 
ing an optically thin correction since this represents a more 
physical upper limit. This gives a value of ^ SCCim „ = 398 
kms -1 (using / = 0.77). The hydrogen Balmer emission is 
likely to be optically thick, as found for the similarly irra- 
diated secondary star in the PCEB NN Ser iParsons et al.l 
l|2010af ). We will show in Section |4~T1 it is the lower limit 
that is the more important. The S/N of the X-shooter data 
is too low to car r y out a similar analysis to that presented 
in lParsons et al.l |2010al ) for NN Ser. A lower limit can be 
obtained by assuming that the emission is uniformly spread 
over the heated face (in reality it is far more likely to be 
brighter towards the substellar point where the heating ef- 
fect is larger). This corresponds to a correction factor of 
/ = 0.42 which gives -/vsec.min = 364 kms -1 . 

Therefore, with the current data we are only able to 
say that the radial velocity semi-amplitude of the centre of 
mass of the secondary star (K s ,, c ) lies somewhere between 
364 kms" 1 and 398 kms" 1 . 



4 DISCUSSION 

4.1 System parameters 

We can combine our light curve model fits, our radial veloc- 
ity measurements and Kepler's third law to determine the 
masses, radii and orbital separation for SDSS J0857+0342. 
These are listed in Table [5] along with the other system pa- 
rameters t hat we have measured , following the procedure 
outlined in IParsons et ail l|2010ah . 

The largest source of uncertainty is K sec since we were 
unable to measure it directly. This translates into a large 
uncertainty in the masses (particularly the white dwarf). 
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Figure 5. Trailed spectrum of the four longest wavelength hydrogen Balmer lines, the scale runs from white (75 per cent of the continuum 
level) to black (120 per cent of the continuum level). Both white dwarf absorption and secondary star emission components are visible 
in each line. 



Table 5. Stellar and binary parameters for SDSS J0857+0342. 
The white dwarf temperature and log g are from the spectroscopic 
fit. zwd is the measured gravitational redshift of the white dwarf. 



Parameter 


Value 


Parameter 


Value 


RA 


08:57:46.18 


Porb (days) 


0.065 096 538(3) 


Dec 


+03:42:55.3 


Twd (K) 


35, 300 ± 400 


u' 


17.731 ±0.011 


WD log g 


7.3 ±0.1 


g' 


17.954 ± 0.006 


Mwd (M ) 


0.514 ±0.049 


r> 


18.256 ±0.007 


M scc (M ) 


0.087 ±0.012 


i' 


18.393 ±0.011 


Rwd (R©) 


0.0247 ± 0.0008 


»' 


18.536 ±0.032 


Rsec (Rq) 


0.1096 ±0.0038 


J 


18.184 ±0.059 


a(R ) 


0.574 ±0.018 


H 


18.618 ±0.160 


i (deg) 


85.5 ±0.2 


K 


18.141 ±0.174 


i\ WD (kms" 1 ) 


64 ±6 


d (pc) 


968 ± 55 


K sac (kms -1 ) 


364 < K < 398 


Sp2 


M8 ± 1 


2wd (kms -1 ) 


18 ±5 



The contours in Figure [6] show the region of allowed mass 
and radius for the white dwarf in SDSS J0857±0342 com- 
pared to other white dwarfs with measured masses and radii. 
Also plotted are mas s-radius relations for carbon -oxygen 
core white dwarfs from lBenvenuto fc Althausl (Il999l ) and he- 
lium core white dwarfs from Panei et al. 1 20071 ) in the same 
temperature range as the white dwarf in SDSS J0857+0342. 
The mass range permitted makes it possible for the white 
dwarf to have either type of core but the measured radius of 
the white dwarf places it far above the carbon-oxygen core 
models, even for very thick hydrogen envelopes. However, it 
is consistent with a helium core provided that the mass is at 
the lower end of the permitted range. Therefore, it is most 
likely that the white dwarf in SDSS J0857±0342 is a helium 
core white dwarf close to the upper mass limit expected from 
evolution (< 0.5Mq). 

The contours in Figure [7J show the region of allowed 
mass and radius for the secondary star in SDSS J0857+0342. 
The mass range places it close to the hydrogen burning limit 
and the radius is consistent with the evolutionary models of 



iBaraffe et al.l (|l99&f ) despite being both heavily irradiated by 
the nearby white dwarf and r apidly rotating. This result is 
consistent with that found bv lParsons et al] (|2010ah for the 
eclipsing PCEB NN Ser, implying that very low-mass stars 
do not appear to be overinflated compared to evolutionary 
models, in contrast to slightly more massive stars (>0.3Mq) 
l|Kraus et alj|20Tlh . The black contours in Figure [7] show the 
region of allowed mass and radius if the white dwarf mass is 
< O.5M~0 (a helium core), in this case the mass and radius 
are slightly smaller (0.08 ± O.O1M and 0.106 ± O.OO2R ). 

The secondary star was detected during the SOFI J 
band eclipse and we measure its magnitude as J scc = 
20.9 ± 0.2 which, at the distance determined for the 
white dwarf, gives an absolute magnitude of Mj iSCC = 
11.0 ± 0.2. This cor responds to a spectral type of M8 ± 1 
l|Hawlev et~aT1 12002ft and is consistent with the measured 
mass (|Baraffe fc Chabrierlll996h . 

Using the X-shooter spectra, the gravitational redshift 
of the white dwarf was found to be 18 ±5 kms -1 . Using the 
measured mass and radius from Table [5] gives a redshift of 
13.2 kms -1 . Correcting the X-shooter measurement for the 
redshift of the secondary star, the difference in transverse 
Doppler shifts and the potential at the secondary star owing 
to the white dwarf gives a value of 12 ± 1 kms -1 which is 
consistent with the measured value. 



4.2 Evolutionary state 

Since SDSS J0857±0342 has emerged from the common en- 
velope stage its evolution has been driven by angular mo- 
mentum loss, which will eventually bring the two stars 
close enough to initiate mass transfer from the secondary 
star to the white dwarf. For short period binaries an- 
gular momentum can be lost v ia gravitational radiation 
l|Kraft et al.lll962l ; lFaulkner|[l97ll ) and a magnetised stellar 
wind, known as magne tic braking (|Verbunt fc Zwaa r] |l98ll ; 
iRappaport et all Il983l ). Magnetic braking becomes much 
weaker in binaries with very low mass stars {M < 0.3Mq) 
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Figur e 6. Mass-ra di us p lot for white dwarfs. Data 
from IProvencal et all fll99Si) . iProvencal et al.l ||2002|) and 
ICasewell et al. I j2009h are plotted. The filled circles are visual 
binaries and the open circles are common proper-motion sys- 
tems. The solid dot s are from the eclipsing PCEBs N N Ser 
llParsons et alj|2010ah and V471 Tau llO'Brien et alj|200ll) . The 
solid black lines correspond to carbon-oxygen core models with 
temperatures of 30,000K and 40,000K with a hydrogen layer 
thickn ess of Mjj/Mwd = 10~ 4 from Bcnvcnut o fc Althauj 
The solid grey lines correspond to helium core mo dels 
with temperatures of 30,000K and 40,000K from [Panei et al.l 
(2007). The dashed line i s the zero-temperature mass -radius 
relation of Eggleton from IVerbunt~fe Rappaport (1988). The 
contours show the 68 and 95 percentile regions of the masses and 
radii for the white dwarf in SDSS J0857+0342. For the measured 
temperature of the white dwarf (35,300K) the measured radius is 
much larger than predicted from the carbon-oxygen core models 
but is consistent with the helium core models. 



since they become fully convective and the magnetic field is 
no longer rooted to the stellar core. However, residual mag- 
netic braking in CVs with fully convective donor stars has 
been suggested both base d on the location of the observed 
orbit al period minimum l|Pattersonl Il998l ; iKolb fc Baraffel 
1999) and on the differences in the temperatures of the 
accretion-heated white dwarfs in strongly ma gnetic and non- 
magnetic CVs (jTownslev fc Gansickel (2009T I. Independent 
evidence for m agnetic braking comes from studies of single 
low mass stars l)Reiners fc Basr3l2008h . 

To determine the past and future evolution of 
SPSS J0857 + 0342 we interpolate the cooling tracks of 
IPanei et al] (|2007r ) for a helium core white dwarf with 
a mass of 0.45 Mq and estimate its cooling age to be 
~ 19.5 x 10 6 yrs, implying that the system has only re- 
cently emerged from the common envelope phase at essen- 
tially the current orbital period. An upper limit on the 
time taken for the system to reach a semi-detached con- 
figuration can be made by assuming that the only mech- 
anism of angular momentum loss is via gravitational radi- 
ation. If this is t he case then mass transfer will begin in 
~ 3.7 x 10 years (jSchreiber fc Gansickei r2003) once the or- 
bital period has been reduced to 66 minutes l|Ritter|[l986T ). 
Since there is likely to be some additional angular mo- 
mentum loss besides gravitational radiation, the true time 
taken to reach a semi-detached configuration is likely to be 
lower. SDSS J0857+0342 will start mass transfer below the 
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Figu re 7. M a ss-rad ius p lot for low mass s tars. Data from 
lLopez-Morales! (|2007}) and iBeattv etldl {2007), also plotted is 
the mass-radius measurement fo r the secondary star in the eclips- 
ing PCEB NN Ser (black point; [Parsons et al.|[2010ah. The solid 
line re presents the theoretical isochrone model from Baraff e~et al.l 
(1998). for an age of 0.4 Gyr, solar metalicity, and mixing length 
a = 1.0, the dashed line is the same but for an age of 4 Gyr. The 
grey contours show the 68 and 95 percentile regions of the masses 
and radii for the secondary star in SDSS J0857+0342. The black 
contours show the same regions if the white dwarf mass is forced 
to be < 0.5M Q . 



observed orbital period minimum for cataclysmic variables 
(CVs), ~ 80 minutes, at which point the white dwarf will 
have cooled to around 13,500K. 

Once the binary becomes mass-transferring the sec- 
ondary star will no longer be in thermal equilibrium and 
the binary will evolve to longer periods and pass into the 
"main" CV population. Therefore, SDSS J0857+0342 is a 
progenitor for a class of CVs yet to be unambiguously ob- 
served: those that began mass transfer from a very low- 
mass star / brown dwarf donor. The few (non-AM CVn) 
CVs observed with periods below the minimum were not 
formed directly from white dwarf / brown dwarf binaries; 
V485 Cen (P orb = 59 minutes) and EI Psc (P or b = 64 
minutes) are both thought to have evolved donor stars in 
which a large fraction of the hy drogen in the core was 
processed prior to m ass transfer dThorstensen et all 120021 : 
iGiinsicke etai1l2003r i. SDSS J1507 (P or b = 66 minutes) was 
though t to have started ma ss transfer with a brown dwarf 
donor (jLittlefair et al.ll2007l ) but recent studies have shown 
that it is more likely to be a halo CV with a lo w-metalicity 
donor |Patterson et al.ll2008l ; IUthas et al]|201ll ). Once it has 
reached a semi-detached state, SDSS J0857+0342 is likely 
to take about a Gyr to evolve back towards a period where 
the s hortest period "main" CVs are found (|Kolb fc Baraffel 
Il999h . 

Our constraints on the system parameters of 
SDSS J0857+0342 make it likely that the progenitor 
for this system lied within, or close to, the brown dwarf 
desert; a paucity of clos e (<5AU) brown dwarf comp anions 
to main-sequence stars jCrether fc Lineweaver|[200q ). 
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5 CONCLUSIONS 

We have used a combination of ULTRACAM and SOFI pho- 
tometry and X-shooter spectroscopy to constrain the sys- 
tem parameters of the 94 minute orbital period eclipsing 
PCEB SDSSJ085746. 18+034255.3. We measure a tempera- 
ture for the white dwarf of 35, 300 ± 400K and a distance of 
968 ± 55pc. By detecting the secondary eclipse in our UL- 
TRACAM photometry we were able to determine an incli- 
nation of i = 85.5° ± 0.2°. We were also able to measure the 
radial velocity amplitude of the white dwarf as 64±6 kms -1 . 
No absorption features from the secondary star were visible 
in the X-shooter spectra, however we were able to determine 
the range of possible K scc values by applying a correction 
factor to the measured radial velocity amplitude of several 
emission lines originating from its heated face. 

By combining the radial velocity measurements with 
the light curve fit we measured the mass and radius of the 
white dwarf to be 0.51 ± O.O5M and 0.0247 ± 0.0008R Q 
respectively. The measured temperature, mass and radius 
are inconsistent with a carbon-oxygen core white dwarf but 
are consistent with a helium core white dwarf at the upper 
mass range for helium core white dwarfs. The secondary star 
has a mass of 0.09 ± O.OIMq and a radius of 0.110 ± 0.004 
consistent with evolutionary models of low-mass stars. Its 
mass places it at the hydrogen burning limit. 

We also find that SDSS J0857+0342 has only recently 
emerged from the common envelope phase and will reach 
a semi-detached configuration in ~ 4 x 10 s years when it 
will become a cataclysmic variable with a 66 minute orbital 
period, at which point its orbital period will increase. 

Due to the detection of the secondary eclipse, 
SDSS J0857+0342 is ideally set up to measure precise masses 
and radii for both of the stars. Currently the limiting fac- 
tor is the poorly constrained K scc . With higher S/N spectra 
a more reliable A'-correction can be made, or potentially 
a direct measurement via absorption lines in the infrared, 
which will vastly improve the precision of the masses and 
allow us to test the mass radius relation for massive helium 
core white dwarfs and very low-mass stars. 
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